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We study multidimensional gravitational models with scalar curvature nonlinearity of the type 
1/R and with form- fields (fluxes) as a matter source. It is assumed that the higher dimensional space- 
time undergoes Freund- Rubin-like spontaneous compactification to a warped product manifold. It is 
shown that for certain parameter regions the model allows for a freezing stabilization of the internal 
space near the positive minimum of the effective potential which plays the role of the positive 
cosmological constant. This cosmological constant provides the observable late-time accelerating 
expansion of the Universe if parameters of the model is fine tuned. Additionally, the effective 
potential has the saddle point. It results in domain walls in the Universe. We show that these 
domain walls do not undergo inflation. 
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I. INTRODUCTION 

There are two great challenges in modern theoret- 
ical physics and cosmology. The first big puzzle con- 
sists in a "dark side" of our Universe. Recent observa- 
tions of the luminosity distances of type la supernovas 
(SNIa), CMB angular temperature fluctuations on de- 
gree scales, and measurements of the power spectrum 
of galaxy clustering indicate (see e.g. [H) that our 
Universe spatially flat with ~ 23% of its critical en- 
ergy in non-relativistic cold dark matter and ~ 73% in 
a smooth component having a large negative pressure 
(dark energy). The latter one results in accelerating 
expansion of our Universe which began approximately 
at redshift z ~ 1 and continues until present time. On 
the other hand, there is also possibility that the late- 
time accelerating expansion of our Universe is caused 
by modification of gravity on Galactic scales. For ex- 
ample, it was proposed H to add a 1/R term in the 
Einstein-Hilbert action to modify General Relativity 1 . 
It is clear that such modification may affect dynamics 
of the Universe at late times of its evolution and on 
large scales where the scalar curvature becomes small. 
In fact, it was shown (see e.g. [J]) that this term can 
provide accelerating expansion of the Universe with- 
out the need of introducing dark energy. 

The second great challenge is possible multidimen- 
sionality of our Universe which naturally follows from 
theories unifying different fundamental interactions 
with gravity, such as string/M-theory. So, there is 
big temptation to explain the dark matter and the 
accelerating expansion of our Universe with the help 
of extra dimensions. However, it is well known that 
dynamical behavior of internal spaces usually results 
in variations of the effective four-dimensional funda- 
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1 Terms with negative powers of curvature can originate due to 
compactification of some fundamental string/M-theory [3|. 



mental "constants" (e.g. gravitational constant, fine 
structure constant, etc.) (see e.g. Q-[£j) and refer- 
ences therein). There are strong experimental bounds 
on such variations [10]. So, one of the main problems 
of higher-dimensional models lies in stable compact- 
ification of the internal spaces. Scale factors of the 
internal spaces play the role of scalar fields moving in 
our four-dimensional space-time. Their dynamics is 
defined by an effective potential in dimcnsionally re- 
duced theory. Thus, the internal spaces are stabilized 
in the case of a minimum of this potential [Til ] . Small 
excitations around this minimum look in our Universe 
as massive scalar fields (gravitational excitons/radions 
[Til ) with Planck scale suppression of their interaction 
with usual matter [7j. Therefore they may play the 
role of dark matter. Additionally, if the minimum 
of the effective potential is positive, it contributes in 
positive cosmological constant providing acceleration 
of the Universe. 

In the present paper, we consider nonlinear gravita- 
tional multidimensional cosmological model with ac- 
tion of the type R + l/R with form-fields as a matter 
source. We also include a bare cosmological term as 
an additional parameter of the theory. It is assumed 
that the corresponding higher-dimensional space-time 
manifold undergos a spontaneous compactification to 
a manifold with warped product structure of the ex- 
ternal and internal spaces. Each of spaces has its 
own scale factor. A model without form-fields and 
bare cosmological constant was considered in paper 
[12] where the internal space freezing stabilization was 
achieved due to negative minimum of the effective po- 
tential. Thus, such model is asymptotically AdS with- 
out accelerating behavior of our Universe. It is well 
known that inclusion of usual matter can uplift po- 
tential to the positive values [lj|. One of the main 
task of our present investigations is to observe such 
uplifting due to the form-fields. Indeed, we demon- 
strate that for certain parameter regions the late-time 
acceleration scenario in our model becomes reachable. 
However, it is not simple uplifting of the negative min- 
imum of the theory [T2] to the positive values. The 
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presence of the form-fields results in much more rich 
structure of the effective potential then in Here, 
we obtain additional branches with extremum points, 
and one of such extremum corresponds to the positive 
minimum of the effective potential. This minimum 
plays the role of the positive cosmological constant. 
With the corresponding fine tuning of the parameters, 
it can provide the late-time accelerating expansion of 
the Universe. Moreover, we show that for this branch 
of the effective potential there is also a saddle point. 
Thus, we obtain domain walls which separate regions 
with different vacua. We demonstrate that these do- 
main walls do not undergo inflation because our ef- 
fective potential is not flat enough around the saddle 
point. 

It is also worth of noting that the effective poten- 
tial in our reduced model has a branchpoint. It gives 
very interesting possibility to investigate transitions 
from one branch to another by analogy with catastro- 
phe theory or similar to phase transitions in statistical 
theory. This idea needs more detail investigation in a 
separate paper. 

The paper is structured as follows. In section [TT] 
we present a brief description of multidimensional 
models with scalar curvature nonlinearity f(R) and 
the form-fields as a matter source. Then, we per- 
form dimensional reduction and obtain effective four- 
dimensional action with effective potential. General 
formulas from this section are applied to our s pecific 
model f(R) — R — fi/R — 2A]j in section UTTl Here, 
we obtain the effective potential minimum cond itions. 
These conditions are analyzed in sections [IV] and |V] 
for the cases of zero and positive effective cosmo logi cal 
constants respectively. Furthermore, in section fVll we 
demonstrate that the positive minimum of the effec- 
tive potential plays the role of the positive cosmolog- 
ical constant and can provide the late-time accelerat- 
ing expansion. Additionally, this minimum is accom- 
panied by a saddle point. It results in non-inflating 
domain walls in the Universe. The main results are 
summarized and discussed in the concluding section 

Em 



II. GENERAL SETUP 

We consider a D — (Do + D')— dimensional nonlin- 
ear gravitational theory with action functional 
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(2.1) 



where f(R) is an arbitrary smooth function of a 
scalar curvature R := R[g] constructed from the 
D— dimensional metric g a \> (a, b — 1,...,D). D' 
is the number of extra dimensions. n 2 D denotes 
the D — dimensional gravitational constant. In action 
(|2.ip , a form field (flux) F has block-orthogonal struc- 
ture consisting of n blocks. Each of these blocks 



is described by its own antisymmetric tensor field 
= 1, . . . , n) of rank di (di-form field strength). 
Additionally, we assume that for the sum of the ranks 
holds £? =1 ^ = D '- 

Following Refs. p2]-[l5[, we can show that the non- 
linear gravitational theory (|2.ip is equivalent to a lin- 
ear theory R — R[g] with conformally transformed 
metric 
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and an additional minimal scalar field 4> = ln[f'(R)]/A 
coupled with fluxes. The scalar field cf> is the result 
and the carrier of the curvature nonlinearity of the 
original theory. Thus, for brevity, we shall refer to the 
field cf> as nonlinearity scalar field. A self-interaction 
potential U(4>) of the scalar field (j> reads 



u(<p) 

where 

A = 



1 



e- B * [ R(ct>)e A * - / 
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Furthermore, we assume that the multidimensional 
space-time manifold undergoes a spontaneous com- 
pactification 



M 



M = Mo X Mi X ... X M„ 



(2.5) 



in accordance with the block-orthogonal structure of 
the field strength F, and that the form fields F«, 
each nested in its own di— dimensional factor space 
Mi (i = 1, . . . , n), respect a generalized Freund- Rubin 
ansatz [Tg| . Here, (Dq — 4)-dimensional space-time 
Mo is treated as our external Universe with metric 

ff (0) (z)- 

This allows us to perform a dimensional reduction of 
our model along the lines of Refs. p| -p3], E3> HI- 
The factor spaces Mj are then Einstein spaces with 

metrics g^ 1 ' = e 2 ^ ( X )^<M which depend only through 
the warp factors a,i(x) :— e 13 '^ on the coordinates 
x of the external space-time Mo- For the corre- 
sponding scalar curvatures holds R py^l = \ l di = T{ 

(in the case of the constant curvature spaces A 1 = 
ki(di — 1), hi — 0, ±1). The warped product of Ein- 
stein spaces leads to a scalar curvature R which de- 
pends only on the coordinate x of the D Q — dimensional 
external space-time Mo: R[g] = R(x). This im- 
plies that the nonlinearity field <f> is also a function 
only of a;: <f> = 4>(x). Additionally, it can be easily 
seen [l2| that the generalized Freund-Rubin ansatz 
results in the following expression for the form-fields: 

(i^)) 2 = f?/a 2d * where p = const. 

In general, the model will allow for several stable 
scale factor configurations (minima in the landscape 
over the space of volume moduli). We choose one 
of them (which we expect to correspond to current 
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evolution stage of our observable Universe), denote 
the corresponding scale factors as and work further 
on with the deviations $ % (x) = fi l (x) — f3 . 

Without loss of generality 2 , we shall consider a 
model with only one di-dimensional internal space. 
After dimensional reduction and subsequent confor- 
mal transformation to the Einstein frame the action 
functional (|2.ip reads 3 
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(2.6) 



where 
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- v / di(D-2)/(D - 2)/3 1 and 
rib/Vd! denotes the -Do - dimensional (four- 
dimensional) gravitational constant. ^exp(di/?o) 
is the volume of the internal space at the present 
time. 

A stable compactification of the internal space Mi 
is ensured when its scale factor cp is frozen at the min- 
imum of the effective potential 



U, 



eff 



--R ie aip + U{<t>) + he cct >e adlip 



, (2-7) 



where R\ := r\ exp (— 2/3q) defines the curvature of the 
internal space at the present time and contribution of 
the form-field into the effective action is described by 
h := K/jj / 2 exp ( — 2(^i/?q) > 0. For brevity we intro- 
duce notations 



a : 



Do 



di 



V di(D-2) ' 

2di - D 
^(D-l)(D-2) 
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III. THE MODEL 

In this section we analyze the conditions of the com- 
pactification for a model with 



f{R)=R- 



It 
R 



2A 



D ■ 



(3.1) 



2 The difference between a generaf modef with n > 1 internal 
spaces and the particular one with n = 1 consists in an addi- 
tional diagonalization of the geometrical moduli excitations. 

3 The equivalency between original higher dimensional and ef- 
fective dimensionally reduced models was investigated in a 
number of papers (see e.g.[l9l|). The origin of this equiva- 
lence results from high symmetry of considered models (i.e. 
because of specific metric ansatz which is defined on the man- 
ifold consisting of direct product of the Einstein spaces). 



Then from the relation f'(R) = exp (A(f>) we obtain 



R = q\ 



? ( e A0 _ 1) 



q = ±1 , s = sign(^) . (3.2) 



Thus, the ranges of variation of <f> are <j> G (-co, 0) for 
H < (s = -1) and 4> G (0, +oo) for fx > (s = +1). 

It is worth of noting that the limit </> — * ±0 (/' — > 
1) corresponds to the transition to a linear theory: 
f(R) — > R — 2A£) and i? — >• R. This is general feature 
of all nonlinear models f(R). For example, in our case 
gU) we obtain f(R) = R(2 - exp (Ac/))) - 2A D — ► 
R — 2A p fo r <j) — ► 0. F rom other hand, for particular 
model (|3.1[) . eq. (|3.2[) shows that the point = 
maps into infinity R, R — ±oo. Thus, in this sense, 
we shall refer to t he p oint = as si ngularity. 

For our model (|3.ip , potential (|2.3[) U (0) reads 
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(3.3) 



It is well known (see e.g. 0,[3]j[i3) that m or ~ 
der to ensure a stabilization and asymptotical freezing 
of the internal space Mi, the effective potential (|2.7p 
should have a minimum with respect to both scalar 
fields ip and 0. We remind that we choose the mini- 
mum position with respect to ip at ip = 0. Addition- 
ally, the eigenvalues of the mass matrix of the coupled 
(<p, 0)-field system, i.e. the Hessian of the effective 
potential at the minimum position, 



J := 



d^U e f f d^Ueff 



(3.4) 



should be positive definite (this condition ensures the 
positiveness of the mass squared of scalar field exci- 
tations). According to the Silvester criterion this is 
equivalent to the condition: 



Jn > , J 22 > , det(J) > 



(3.5) 



It is convenient in further consideration to introduce 
the following notations: 



0o ■= (f>\extr , X := V se A <t>° - s > -> A (s= _ 1} < 1 . 

(3.6) 

Then we can rewrite potentials U((f>), U e f /(</?, 0) and 
derivatives of the U e f / at an extremum (possible min- 
imum) position (tp = O,0o) as follows: 

u = u\ extr = (i + sX 2 )~ B/A (qsy/\}T\X + a d) , 



Uefflextr = —Rl + U (X) + h(l + sX 2 ) C/A 



(3.7) 
, (3.8) 



dJJ t 



<P e ff \Gxtr 



-Ri + bU (X) 



+ (adi + b)h (1 + sX 2 ) c/A = 0, (3.9) 
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d^U eff \ extr = ch(l + sX 2 ) c/ -BUq(X) 

+ ^(l + ^ 2 r SM = 0,(3.10) 

d 2 w U eff \ extr = - { ^±^R 1+ b 2 U (X) 

+ {ad 1 +b) 2 h(l + sX 2 ) c/A , (3.11) 



to a minimum of the effective potential. With respect 
to this solution we define parameters in the relation 
(|4.1[) . Therefore, we must distinguish now which of X p 
corresponds to the minimum of U e ff. Let us investi- 
gate solutions (|4.3[) for the purpose of their satisfac- 
tions to conditions e^ > 0, U (X) > and X p > 0. 

The condition e A ^° = 1 + sX 2 > 0: 

Simple analysis shows that solutions X p satisfy this 
inequality for the following combinations of parame- 
ters: 



d y<tPeff\extr = chad 1 (l + sX' 



,c/A 



(3.12) 



eff\extr = ch (c - A + 2B) (l + sX 2 , 

+ B(A - B)Uo(X) 

qSy^A 2 , 2 s2-B/A 

AX^ {1 + SX > 



c/A 



. (3.13) 



The most natural strategy for extracting detailed in- 
formation about the location of stability region of pa- 
rameters in which compactification is possible would 
consist in solving (|3.10[) for X with subsequent back- 
substitution of the found roots into the inequalities 
(|3.5p and the equation (|3.9[) . To get the main features 
of the model under consideration, it is sufficient to in- 
vestigate two particular nontrivial situations. Both of 
these cases are easy to handle analytically. 



n > 

(s = +1) 


p = ±1 : z e (-oo.+oo) 


At < 

(s = -1) 


p = +l : z e (-2, -z Q ) U (z ,+oo) 
p = -l : z e (-00,-20) U (zo,2) 

Z = Zq 



(4.4) 

The condition U (X) > 0: 

As appears from eq. (|3 7|) , this condition takes place 
if X p satisfies inequality 2qsX p + z > which leads to 
the conditions: 



n> 
(s = +1) 


p = +1 : z G (— 00, +00) 


n < 
(s = -1) 


p = +1 : z e (z , +00) 
p= -1 : z£ (-20,2] 
z = z 



(4.5) 



The condition A p > 0: 

This condition is satisfied for the combinations: 



IV. ZERO EFFECTIVE COSMOLOGICAL 
CONSTANT: A eff = 

It can be easily seen from eqs. (|3.8j) and (|3.9p that 
condition A e ff = f/ e //|extr = results in relations 

iZi = 2d 1 h (1 + s* 2 ) cM = [/•„(*) . di > 2 , 

(4.1) 

which enable us to get from eq. (13.10P quadratic equa- 
tion for X 

(di + l)A 2 +qsdizX - s(di-l) = 0, z = 2Az>/V4VI 

(4.2) 

with the following solutions: 

Xp = ^2^1) (~ z+p \f 2+4 4^) ' 

p = ±1. (4.3) 

In the case s = — 1 parameter z should satisfy con- 
dition |z| > z = 2^/d 2 — 1/rfi < 2 and for z = z 
two solutions X p degenerate into one: X p = Xo = 

-qs^{d 1 - l)/{d 1 + 1). 

Because of con ditio ns /i > and e A< ^° = 1 + sX 2 > 
0, the relations (|4.1[) show that parameters i?i and 
£/oPO should be non-negative: Ri > 0, E/oPO > 0. 
Obviously, only one of the solutions (|4.3p corresponds 



M > 

(» = +l) 


<? = +l : < 
9 =-l:< 


' z < 
z > 

' z < 
[ z > 


A + > A_ < 
X+ > x_ < 

A+ < A_ > 
A+ < X_ > 


At < 

(,s = -1) 


9 = +l : < 
q=-l:< 


1 z < 
z > 

> z < 
z > 


X+ < A_ < 
A+ > A_ > 

A+ > A_ > 
X+ < x_ < 



(4.6) 

The comparison of (|4.4p . (|4~5)) and l|4.6p shows that 
they are simultaneously satisfied only for the following 



combinations: 



At > 
0» - +1) 


A3 = +1 : ? = +l : z€ (—00, +00) 


At < 
(*=-!) 


p=+l : 9 = +1 : z <E (z , +00) 
p=-l : q = +1 : z G (z ,2) 
g = +1 : z = zq 



(4.7) 



Additionally, the extremum solutions X p should 
correspond to the minimum of U e f / ■ The inequalities 
(|3.5p are the sufficient and necessary conditions for 
that. We analyze them in the case of four-dimensional 
external space D Q — 4. Taking into account defini- 
tions ((231), (J2ll), (|33T]) - (f3~T3l) and relations lpL"T|) . for 
J11, J22 and J21 we get respectively: 

Jn = J^MX P ) , (4.8) 
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6di(2qsX p 
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dl+4 

sXT] d i+ 2 



(di - l)(di + 2)(di + 3) 
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4X|(di + 3) 



(1 + SX2)W 



(4.9) 



J 2 i = - 



(di - A)(2qsX p + z) 
(di-l)(di+2) 



di(di +3) 



dl+4 

sX p 2 ) dl+2 



. (4.10) 



We supposed in these equations that each of X p can 
define zero minimum of U e ff- In what follows, we 
shall check this assumption for every X p with corre- 
sponding combinations of signs of the parameters s 
and q in accordance with the table (|4.7[) . 

According to the Silves ter c riterio n (13. 5|) . Jn should 
be positive. Thus eqs. (|4.1|) and (|4.8p result in the 
following conclusions: the potential Uo should be pos- 
itive Uq > 0, the internal space should have positive 
curvature R\ > (hence, d± > 1) and its stabilization 
(with zero minimum A e ff = 0) takes place only in 
the present of form-field (h > 0). Transition from the 
non-negativity condition Uo > to the positivity one 
Uo > corresponds to the only substitution in (|4.5p 
(z ,2] — > (zo,2) for the case s — —1, p = —1. Ex- 
actly this interval (zq, 2) appears in concluding table 
(|4.7p . Therefore, Jn is positive for all X p from the 
table (|4~7) . 

Concerning expressions J22 and det(J) = Jn J22 — 
J^2, graphical plotting (see Fig [T] and FigJ5]) demon- 
strates that they are negative for s — +1, p = 
+1, g = +1 and s = —1, p = —1, g = +1 but positive 
in the case s = —1, p = +1, q = +1. For this latter 
combination 2 € (zo,oo)- The case s = —1, g = +1 
and z = zo should be investigated separately. Here, 

X p = Xo = \/(di 
we obtain: 



l)/(di + 1) and for J22 and J21 



H\X- 3 (1 - 
12(di - 1) 



sXl) 



" d t +2 



(di + 2) 2 



J21 



(d 1 + l) 2 (d 1 +2)(di+3) 



n\x^{i + sx%y d ^ 

2\/2(di -4)(di - 1) 



(4.11) 



(di + l) 2 (di+2) v /di(d 1 +3) 



(4.12) 



It can be easily seen from eqs. (|4.11j) and (|4.12p that 
J 22 > for di ^ 4 and J 22 = J21 = for di = 4. 
Additionally, det(J) > for d x ^ 4. 

Thus, we can finally conclude that zero minimum of 
the effective potential t/ e // takes place either for ,s = 
— 1, q = +1, z G (zo,oo),Vdi > 1 (position of this 
minimum is defined by solution (|4.3p with p = +1) or 



for s = —1, q = +1, z = zo,Vdi ^ 4. Concerning 
the signs of parameters, we obtain that [i < and 
A D > 0. 



f I s=+l p=+l q=+l 




FIG. 1: Typical form of J22/VH ( e< l- EU) f° r param- 
eters s = +1, q = +1, p = +1 (upper) and s = —1, q 
■ I . p = —1 (lower). For both of these combinations of 
the parameters, J22 < 0. 



V. DECOUPLING OF EXCITATIONS: di = D 

It can be easily seen from eq. (|2.8p that in the 
case di = Z?o parameter c = that leads to condition 
dtd,Ueff\extr = (see eq. (|3.12p . Thus the Hessian 
(|3.4p is diagonalized. It means that the excitations 
of the fields if and <f> near the extremum position are 
decoupled from each other 4 . 

Dropping the h term in eq. (|3.10p (because of c = 0) 
and taking into account eq. (|3 7p , we obtain quadratic 
equation for X 



(D + 2)X 2 + DqszX - (D - 2)s = , 



(5.1) 



which for Do — d\ exactly coincides with eq. (14. 2p . 
Thus, in spite of the fact that we does not use the 



4 In the vicinity of a minimum of the effective potential, 
squared masses of these excitations are = Jn and 

ml = J 22 . 
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FIG. 2: Typical form of J22/VM ( u PP er ) an d det(J)/|/i| 
(lower) (see eqs. (|4.8[) - (14.10[0 for parameters s — — 1, g = 
+ 1, p = +1. For this combination of the parameters, 
J22 > and det J > 0. 



condition A e // = 0, we obtain in the case d\ = Dq 
precisely the same solutions (|4.3j) . However, parame- 
ters i?i , Uq and h satisfy now relations different from 
(14. ip . For example, for the most physically interesting 
case Dq = d\ = 4, eqs. (|3.0|) and (|3.8|) result in the 
following relations: 



Ri 



U (X) 



A eff (X) = ±U (X)-h. 



Nonzero components of the Hessian read 
J11 = l[9h-U Q (X)} , 



(5.2) 



(5.3) 



'22 



4(z + 2qsX p ) 3qs{l + sX p 2 ) 2/3 



21(1 + sX2)4/3 



14X3 



In this section we are looking for a positive minimum 
of the effective potential. It means that A e // > 0, 
Jn > and J22 > 0. From the positivity of Jn and 
A e ff we obtain respectively 5 : 

J u >0: l6h>Ri >WU (X)/9>8A eff (5.4) 



and 

Ae// > 



h > fix/16 > U a (X)/9 >h/3>0 



These inequalities show that for the considered model 
positive minimum of the effective potential is possible 
only in the case of positive curvature of the internal 
space i?i > and in the presence of the form field 
(h > 0). 

To realize which combination of parameters s,p and 
q ensures the minimum of the effective potential, we 
should perform analysis as in the previous case with 
Ae// = 0. However, there is no need to perform such 
analysis here because solutions of eq. (|5.ip coincides 
with (|4.3[) and all conditions for X p and C/ are the 
same as in the previous section. Thus, we obtain con- 
cluding table of the form (|4.T[) . Additio nally , it ca n be 
easily seen that expressions of J22 in (|5.3[) and (14. 9p 
exactly coincide with each other if we put d\ = 4 in 
the latter equation 6 . Hence, we can use FigfT] and 
FigH] (for the lines with d\ — 4) to analyze the sign of 
J22 ■ With the help of these pictures as well as keeping 
in the mind that J2i(d\ = 4, z = zq) = 0, we obtain 
that the only combination which ensures the positive 
minimum of U e ff is: s = — 1 , p = +1 , q = +1 and 
z G (zq = 1.9 3 6, + 00). It is clear that potential U q(X) 
in eq s. (|5.2p - (|5.5p is defined by solution of eq. (|5.1[) 
(i.e. (|4.3p for d\ = 4) with this combination of the pa- 
rameters. Because s = — 1 and z > , the parameters 
fi and Ad should have the following signs: /1 < and 
A D > 0. 

Additionally, it is easy to verify that second solution 
of eq. (|5~Tj) X_ (with p = -1 and s = -1, q = +1) 
does not correspond to the maximum of the effective 
potential U e ff. Indeed, we have here d^,U e ff(X^) = 

but d v Ueff{X_) ^0. 

Fig. [5] demonstrates the typical profile ip = of 
the effective potential U e ff((p,4>) in the case of posi- 
tive minimum of t/ e // considered in the present sec- 
tion. This picture is in good concordance with the 
table (|4.7p . According to this table, positive extrema 
of U e ff are possible only for the branch q = +1 of 
the solution (|3.2p (solid lines in Fig. [3]). We see that 
for Zq < z < 2 we can have 3 such extrema: one 
for positive fj, > — > s = +1 and two for negative 
fj, < — >s = — 1. Our investigations show that in 
the left half plane (i.e. for /1 < 0) the right extremum 
(p = +1 in eq. (|4.3p ) is the local minimum 7 and the 



It is interesting to note that (in the case Dq = d\ = 4) rela- 



tions l|5.1j l. Jn in l|5.3|l and inequalities l|5.4| l. I|5.5j l coincide 
with the analogous expressions in paper .13] with quadratic 
nonlinear model. This is not surprising because they do not 
depend on the form of nonlinearity f(R) (and, consequently, 
on the form of {7(0)). However, the expressions for J22 are 
different because here we use the exact form of U (<j>) . 

6 It follows from the fact that in l|4,9j l we already put Do = 4. 
Although we use in this equation the relation 114.11 1 between 
h and Uo, it enters here in the combination which is propor- 
tional to c. Thus, this combination does not contribute if we 
put additionally di = 4. 

7 It can be easily seen for the branch q = +1, s = — 1 that 
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left maximum (p = — 1) is not the extremum of U e ff 
because here d v U e ff 7^ 0. Analogously, maximum 
in the right half plane /i > (which corresponds to 
(p = +l)-solution Q4.3p ) is not the extremum of U e ff. 
For completeness of picture, we also included lines cor- 
responding to the branch q = — 1 (dashed lines in Fig. 
[3]). The minimum of the right dashed line (for X p 
with s = +1, p = —1, q = —1) does not describe the 
extremum of U e ff because again d v U e ff ^ 0. 



fi<o u eff ,u>o 



\ 




/ -4 -2 


\ \ 2 4 


/ -0.5 




-1 





FIG. 3: Profile ip = of the effective potential U e ff(ip, <f>) 
for parameters z = 1.99 (A D = 1.99/2 , = 1), D = 
di = 4 and h — Up /6. The rest of parameters can be 
found from relations (|5.2[1 . Solid and dashed lines describe 
branches q — +1 and q = — 1 respectively. For the solid 
line, there is only one local minimum of U e j[ which is 
defined by solution X p with the following parameters: s — 
— 1, p = +1, q = +1. Left (s = -1, p = -1, q = 
+1) and right (s = +1, p = +1, q = +1) maxima are 
not extrema of U e ff because in these points d$U e ff = 
but d v U e ff 7^ 0. Analogously, for right dashed line the 
solution X p (with s = + 1, p = —1, q = —1) does not 
correspond to the extremum of U e ff. 



VI. COSMIC ACCELERATION AND 
DOMAIN WALLS 

Let us consider again the model with d\ = Dq = 4 
in order to define stages of the accelerating expansion 
of our Universe. It was proven that for certain con- 
ditions (see (|5.2p - (l5.5[) ) the effective potential Ueff 
has local (for z < z < 2) or global (for z > 2) 
positive minimum. The position of this minimum is 



U($ — > — oo) — » +oo for z > 2 and U(<j> —* — oo) — > — oo for 
z < 2. Thus, for z > 2 this minimum becomes global one. 



(ip = 0, <)> = {l /A) ln(l + sX p 2 )) where s = -1 , p = 
+ 1 , g = +1 and z G (z ,+oo). Obviously, posi- 
tive minimum of the effective potential plays the role 
of the positive cosmological constant. Therefore, the 
Universe undergoes the accelerating expansion in this 
position. Thus, we can "kill two birds with one stone" : 
to achieve the stable compactification of the internal 
space and to get the accelerating expansion of our ex- 
ternal space. 

We associate this acceleration with the late-time ac- 
celerating expansion of our Universe. As it follows 
from cqs. (|5.2p and (|5.5|) . positive minimum takes 
place if the parameters are positive and the same order 
of magnitude: A e ff ~ Ri ~ U(X) ~ h > 0. On the 
other hand, in KK models the size of extra dimensions 
at present time should be 6(o)i < lCU 17 cm ~ ITcV -1 . 
In this case Ri > — 10 34 cm~ 2 . Thus, for the 

TeV scale of 6( )i ~ ITeV we get that A e ff ~ Ri ~ 
U(X) ~ h ~ ITeV 2 . Moreover, in the case of natural 

condition ~ v|m| we obtain that the masses of 
excitations m v ~ ~ ITeV. The above estimates 
clearly demonstrate the typical problem of the sta- 
ble compactification in multidimensional cosmological 
models because for the effective cosmological constant 
we obtain a value which is in many orders of magni- 
tude greater than observable at the present time dark 
energy ~ 10 -57 cm~ 2 . The necessary small value of 
the effective cosmological constant can be achieved 
only if the parameters R\ , U(X) , h are extremely fine 
tuned with each other to provide the observed small 
value from equation A e ff(X) = Uq(X)/3 — h. We 
see two possibilities to avoid this problem. Firstly, 
the inclusion of different form-fields / fluxes may result 
in a big number of minima (landscape) [20]- 23] with 
sufficient large probability to find oneself in a dark en- 
ergy minimum. Secondly, we can avoid the restriction 
Ri ~ b7o)i ~ 10 34 cm~ 2 if the internal space is Ricci- 
flat: Ri = 0. For example, the internal factor-space 
Mi can be an orbifold with branes in fixed points (see 
corresponding discussion in [24|). 

The WMAP three year data as well as CMB data 
are consistent with wide range of possible inflation- 
ary models (see e.g. Therefore, it is of interest 
to get the stage of early inflation in our model. It 
is well known that it is rather difficult to construct 
inflationary models from multidimensional cosmolog- 
ical models and string theories. The main reason of 
it consists in the form of the effective potential which 
is a combination of exponential functions (see e.g. eq. 
Q2.7p ). Usually, degrees of these exponents are too 
large to result in sufficiently small slow-roll parame- 
ters (see e.g. O)- Nevertheless, there is a possibility 
that in the vicinity of maximum or saddle points the 
effective potential is flat enough to produce the topo- 
logical inflation [25|, [26|, [27j • Let us investigate this 
possibility for our model. 

As stated above, the value if = corresponds to 
the internal space value at the present time. Following 
this statement, we found the minimum of the effective 
potential at this value of <p. Obviously, the effective 
potential can also have extrema at ip ^ 0. Let us 
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investigate this possibility for the model with d\ = 
Dq = 4, i.e. for c = 0. In this case, the extremum 
condition of the effective potential reads 



1 



b)e {a+b)vo + bU e btpo 
{ad x + b)he {adl+b ^° = (6.1) 



and 



d s U t 



<t> U eff\ ip0y< f >( 



= e bvo 



dU 



= 



= 



00 



dU 

0o > °* 

(6-2) 

Here, ipo and c/)q define the extremum position. It 
clearly follows from eq. (|6.2|) that <fio is defined by 
equation dll/d(f> = which does not depend on ip. 
Therefore, extrema of the effective potential may take 
place only for 0o which correspond to the solutions 
X p of eq. (|5.1|) and different possible extrema should 
lie on the sections X p — const. So, we take X + (with 
s = — 1 and q = +1) which defines the minimum of 
U e ff in the previous section. Hence, Uo in eq. (|6.ip 
is the same as for eq. (I5.2p . 

Let us define no w ipp from eq. (|6.1[) . With the 
help of inequalities f|5.4|) and (|5.5p we can write ft = 
nf/o wher e n € (1/ 9, 1/ 3). Taking also into account 
relations (|5.2[) . eq. (|6.1[) can be written as 



V 



1 

3?; 



1 

3/? 



0. 



(6.3) 



where we introduced the definition y = exp(apo) = 

exp (tpo /v3) and p ut d \ = Dq = 4. Because y = 1 is 
the solution of eq. (|6.3[) . the remaining three solutions 
satisfy the following cubic equation: 



y 3 + y 2 + y - — = . 

in 



(6.4) 



It can be easily verified that the only real solution of 
this equation is 



where 



1 ( 1 

2/o = r -1-2I/ + - 
3 \ ^ 



2 l/3, 



(6.5) 



(9n 2 + 7n 3 + 3\/9n 4 + Un 6 + 9n 6 ) 
1 1 
9' 3 



1/3 



(6.6) 



Thus (fio(yo) and 4>q(X + ) define new extremum of 
U e f f . To clarify the type of this extremum we should 
check signs of the second derivatives of the effective 
potential in this point. First of all we should remem- 
ber that in the case c = mixed second derivative dis- 
appears. Concerning second derivative with respect to 
d>, we obtain 



J: 



d 2 U, 



22 



eff 



= e bvo 



d 2 u 



> 



(6.7) 



because we got in previous section d 2 U / d<fi 2 
0. Second derivative with respect to ip reads 



00 (-^4 



> 



Jll 



d 2 U eff 



dip 2 

12ny4] = 2y 2 U 



Vo,<Po 

,1.1 - 0„,27 



y 2 u 

(3n + l)y 



(6.8) 



where we took into account eq. (|6.3p . Simple analysis 
shows that (3n+ l)y < 4/3 for n & (1/9, 1/3). Keep- 
ing in mind that Uq > we obtain Jn < 0. Therefore, 
our extremum is the saddle surface 8 . Figure |4] demon- 
strates contour plot of the effective potential in the 
vicinity of the local minimum and the saddle point. 




FIG. 4: Contour plot of the effective potential U e ff{y>,4>) 
for parameters z = 1.99 (Ad = 1.99/2 , |^| = 1), D = 
di — 4 an d h — Uo/6. The rest of parameters follows from 
relations (|5.2p . We choose the branch corresponding to 
s — — 1, q = +1. This plot clearly shows the minimum 
and the saddle points of the effective potential. 

Therefore, we arrived at very interesting possibility 
for the production of an inflating domain wall in the 
vicinity of the saddle point. The mechanism for the 
production of the domain walls is the following [H| . If 
the scalar field ip is randomly distributed, some part 
of the Universe will roll down to ip — 0, while in oth- 
ers parts it will run away to infinity. Between any two 
such regions there will appear domain walls. In Ref. 
[26], it was shown for the case of a double- well po- 
tential Vdw{<p) = (A/4)(</? 2 — £ 2 ) 2 that a domain wall 
will undergo inflation if the distance £ between the 
minimum and the maximum of Vdw exceeds a criti- 
cal value £ cr = 0.33Afp; — * Ko^ cr = 1.65. In our 



00 



8 Similar analysis performed for the branch with s = +1, q = 
— 1 (right dashed line in the Fig. |3j shows the existence of 
the global negative minimum with ip ^ along the section 
X— = const. 



case it means that the distance \(fo\ between the lo- 
cal minimum and the saddle point should be greater 
than £ cr : \cpo\ > 1.65. Unfortunately, for our model 
(po(n) < 1.65 if n S (1/9,1.3). For example, in the 
most interesting case n — > 1/3 (where A e // — * (see 
eq. (!5.2p D we obtain \(p$\ — > 1.055 which is less than 
£ cr . Moreover, our domain wall is not thick enough 
in comparison with the Hubble radius. The ration of 
the characteristic thickness of the wall to the horizon 
scale is given by r w H ps \U eff /3d vv tf e //|^ (x+ ) -> 
0.454 for n — ■> 1/3 which is less than the critical 
value 0.48 for a double-well potential. Thus, here 
there is no a sufficiently large (for inflation) quasi- 
homogeneous region of the energy density. Our po- 
tential is too steep. Obviously, the slaw roll param- 
eter e fa (1/2) (d v U eff /U eff )l o rt>o is equal to zero in 
the saddle point. However, another slow roll parame- 
ter rj « |c^£/ eJ/ /£/ e// |^ o -> -1.617 for n - 1/3. 

Therefore, our domain walls do not inflate in contrast 
to the case R 4 in Ref. [27|]. 

In Fig. [5] we present comparison between our po- 
tential (solid line) and a double- well potential (dashed 
line) in the case n = (1 — 0.001)/3. We see that our 
potential is flatter than a double- well potential around 
the saddle point. However, our calculations show that 
it is not enough for inflation. 
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FIG. 5: Comparison of the potential U e f f(tp,<j>a) with a 
double well potential for parameters z = 1.99 (Ad = 
1.99/2 , \p\ = 1), Do = di = 4 and n = (1 - 0.001)/3. 



VII. CONCLUSIONS AND DISCUSSIONS 

We have shown that positive minimum of the ef- 
fective potential plays the double role in our model. 
Firstly, it provides the freezing stabilization of the 
internal spaces which enables to avoid the problem 
of the fundamental constant variation in multidi- 
mensional models (@-[l|). Secondly, it ensures the 
stage of the cosmic acceleration. However, to get the 
present-day accelerating expansion, the parameters of 
the model should be fine tuned. Maybe, this problem 
can be resolved with the help of the idea of landscape 



of vacua ( 20]-[23]). We intend to investigate this pos- 
sibility in our forthcoming paper. 

We have additionally found that our effective poten- 
tial has the saddle point. It results in domain walls 
which separates regions with different vacua in the 
Universe. These domain walls do not undergo infla- 
tion because the effective potential is not flat enough 
around the saddle point. 

It is worth of noting that minimum in Fig. 3 (left 
solid line) is metastable. In other words, classically it 
is stable but there is a possibility for quantum tun- 
nelling both in <j) and in (p directions (see Fig. 4). 
We can avoid this problem in tj> direction in the case 
of parameters z > 2 (see footnote 7). However, tun- 
nelling in ip direction (through the saddle) is still valid 
because U e ff(<p,(f>o) pa e bv U((j)o) — > for ip — * — oo 
which is less than any positive A e f f . It may result in 
the materialization of bubbles of the new phase in the 
metastable one (see e.g. [13 )■ Thus, late-time accel- 
eration is possible only if characteristic lifetime of the 
metastable stage is greater than the age of the Uni- 
verse. Careful investigation of this problem (includ- 
ing gravitational effects) is rather laborious task which 
needs a separate consideration. As we mentioned in 
footnote 8, there is also the global negative minimum 
for right dashed line in Fig. 3 (it corresponds to the 
point (ip = 0.67, (j> = 1-66) for parameters taken in 
Fig. 3). This minimum is stable both in classical and 
quantum limits. However, the acceleration is absent 
because of its negativity. 

Another very interesting feature of the model un- 
der consideration consists in multi- valued form of the 
effective potential. As it can be easily seen from eqs. 
(|2.7p and (|3 ,3|) , for each choice of parameter /i poten- 
tial U(4>) (and consequently U e ff) has two branches 
(q = ±1) which joint smoothly with each other at 
= (see Fig. [3|). It gives very interesting possibility 
to investigate transitions from one branch to another 
one by analogy with catastrophe theory or similar to 
the phase transitions in statistical theory. However, as 
we mentioned above, in our particular model the point 
= corresponds to the singularity R,R — > ±oo. 
Thus, the analog of the second order smooth phase 
transition through the point cj> = is impossible in 
our model. Nevertheless, there is still a possibility for 
the analog of the first order transition via quantum 
jumps from one branch to another one. In what fol- 
lows, we plan to investigate such "phase transitions" 
for non-linear multidimensional models f{R). 

To complete the paper, we investigate some limit- 
ing cases. Firstly, we consider the limit h — > (for 
arbitrary Dq and d\) where the form-fields are ab- 
sent. From eqs. (|3.7I) - (|3.13[) we obtain the following 
system of equations: 

Ri = -^-rU (X), U e ff\ extr — -——iU (X) (7.1) 



a + b 



and 

J n = -abU {X) 7 J 21 =0, 
sAB 



J22 — 



B{A-B)-^{l + sX') 



U (X). (7.2) 
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Since for minimum should hold true the condition 
Jn > 0, we arrive at the conclusion: R\, Uo, U e ff < 
0. Consequently, the minimum of the effective poten- 
tial as well as the effective cosmological constant is 
negative and accelerating expansion is absent in this 
limit. Therefore, the presence of the form-fields is the 
necessary condition for the acceleration of the Uni- 
verse in the position of the freezing stabilization of 
the internal spaces. Additionally, it can be easily seen 
that the extremum position equation takes the same 
form as (|5.ip . Simple analysis show that minimum 



takes place for the branch: s = +1 (i.e. /i > 0), 
p = — 1, q = —1 and z S (— oo, +oo). If additionally 
we demand z — > (i.e. Ad — ► and /i is fixed) then 
we reproduce the results of Ref. [l2| ■ 
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